arXiv:1503.01010v3 [quant-ph] 11 Sep 2015 


Quantum simulations of dissipative dynamics: time-dependence instead of size 
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The simulation of quantum systems has been a key aim of quantum technologies for decades, and 
the generalisation to open systems is necessary to include physically realistic systems. We introduce 
an approach for quantum simulations of open system dynamics in terms of an environment of mini¬ 
mal size and a time-dependent Hamiltonian. This enables the implementation of a continuous-time 
simulation with a finite environment, whereas state of the art methods require an infinite environ¬ 
ment or only match the simulation at discrete times. We find necessary and sufficient conditions 
for this Hamiltonian to be well behaved and, when these are not met, we show that there exists an 
approximate Hamiltonian that is, and look into its applications. 


Introduction 

Every quantum system inevitably interacts with its en¬ 
vironment. As quantum simulations are a key aim of 
quantum technologies [1] , the question of how open sys¬ 
tems can be simulated efficiently on a quantum computer 
is one which has received significant research interest in 
recent years W- In addition to this, open systems have 
been shown to be useful for state engineering p HU], 
and as an alternative model of quantum computation 
m- The dynamics of an open quantum system is de¬ 
scribed by a time-dependent quantum channel, which is 
a completely-positive trace-preserving map which acts on 
quantum states in such a way that probabilities stay well 
defined [MIS]- Such channels are often derived by as¬ 
suming that the system interacts via a Hamiltonian with 
an environment, which is then traced out because it is not 
experimentally accessible. A possible way of simulating 
open systems is therefore to recreate this system-bath in¬ 
teraction in a controllable manner, called a dilation, and 
implement the dynamics directly [16] . However, this is 
rarely feasible due to the infinite size of the environment 
and the intrinsic difficulty in engineering such a system. 
Another method is to use Stinespring’s theorem [13 EH], 
which states that every quantum channel at a given point 
in time is equivalent to a unitary acting on a larger state 
followed by tracing out the ancilla, to create a finite di¬ 
lation. This, however, has the disadvantage of modelling 
the evolution to a fixed point in time only, rather than 
replicating the dynamics for all times. Performing a se¬ 
ries of Stinespring dilations one after the other allows the 
evolution to be matched at discrete times, at the cost of 
a large increase in the ancilla space Pi- 

The central idea of this paper is to find, for a time- 
dependent quantum channel which describes the evolu¬ 
tion of a system, a finite dilation such that the evolution 
matches at all times, and is smooth so that it gives rise 
to a well behaved Hamiltonian. Having such a dilation 
allows the dynamics of an arbitrary open system to be 
simulated continuously in time simply by acting on a fi¬ 


nite system with a physically sensible Hamiltonian. This 
is useful for quantum simulations, particularly when the 
time at which the system will be measured is not known 
beforehand, such as in schemes which rely on photon 
counting. It is also applicable in cases where we desire to 
monitor the system continuously via weak measurements 
m- This allows information about the behaviour of the 
system over an interval of time to be recovered; a situa¬ 
tion where prior approaches which evolve the system to 
a fixed point in time would fail. Furthermore, by shifting 
the complexity from inhnite space to time-dependence, it 
provides a model on which open systems can be studied 
easily [UU], which we use to investigate how the system- 
environment interaction is affected by adding controls to 
the system. 

It is known that given two channels which are close 
by it is possible to find two unitary dilations which are 
also close by jUTj; we add the stronger constraint of 
the unitaries varying smoothly so that we can define a 
Hamiltonian for the dilation. We provide an explicit 
method for constructing such dilations, and establish pre¬ 
cise relations between the continuity and boundedness 
of this Hamiltonian and properties of the original quan¬ 
tum channel. We show that it is always possible to find 
such a Hamiltonian which matches the dynamics arbi¬ 
trarily well, provided that the original evolution is an¬ 
alytic in time. Although our methods are not limited 
to Markovian channels, we pay particular attention to 
these as they constitute some of the most common types 
of noise encountered in quantum information [22] . These 
are the class of channels which are memoryless so that 
they can be expressed in terms of an equation of motion, 
p = Lt (p ), where Lt is a Lindbladian [Til I23H25] . This 
allows us to understand this work as raising dilations 
to the levels of generators; mapping Lindbladians into 
Hamiltonians. More generally, we study the usefulness 
of these dilations in the simulation of open systems, give 
several explicit examples for common quantum channels, 
and look at the effect of adding Hamiltonian controls; as 
well as how to generalise the approach to systems con- 
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sisting of many qubits. 


Method 

We begin by presenting a series of steps, illustrated in 
FIG.[^ which allows us to construct a dilation. Although 
the Stinespring dilation is highly non-unique, we use this 
method for two reasons. Firstly, it gives us an explicit 
way of constructing a unitary dilation [UlinillZl- Sec¬ 
ondly, and for our purposes crucially, it enables us to 
study the properties of the Hamiltonian by letting us fol¬ 
low the time-dependency of the objects throughout their 
transformations. 



FIG. 1. A schematic of how to construct the dilation from 
a quantum channel. All the steps preserve the analyticity in 
time apart from separating the diagonalised Choi state into 
Kraus operators which introduces a square root. It is differ¬ 
entiating this square root that introduces the possibility of 
discontinuities and divergences. 

We start with a family of quantum channels, et('), 
which is analytic in time. In the case of Markovian chan¬ 
nels, these are generated by a Lindbladian according to 
£((•) = T exp ^Jp clrLT^-)^ where T is the time-ordering 
operator, but we do not limit ourselves to such cases. 
The next step is to construct the (unnormalised) Choi 
state which arises from the Choi-Jamiolkowski isomor¬ 
phism [18], and is given by A(t) = {ct ®X) jH) (fl| where 
|r2) = \jj) is a maximally entangled state between 
the system’s Hilbert space and its duplicate, and I is the 
identity map. This is equivalent to reshuffling the ele¬ 
ments of the channel represented as a matrix |28j . This 
is just a linear transformation, therefore A(t) is also an¬ 
alytic in t. As A(t) is Hermitian and positive, it can be 
decomposed into its eigenvalues and eigenvectors, Xk{t) 
and \vk(t)) respectively, where the index k runs from 
1 to the Kraus rank of the channel, R, which is upper 
bounded by cP. It is known via perturbation theory 
that, as A(t) is Hermitian and analytic, its eigenvalues 
and vectors are also analytic for real t [5^. This al¬ 
lows us to write the channel in its Kraus representation, 


£t(p) = ^k{t) pMl{t), where the Kraus operators are 

Mkit) = ^ sJXkit) {i ® Xl\vk{t)®j) \i) (j| (1) 

b 

for any choice of basis {|i)}. The presence of the square 
root here is key. It results in the Kraus operators be¬ 
ing continuous everywhere, and smooth everywhere apart 
possibly from individual points. 

The final step in finding a dilation is solving 
X) = TYB[U{t) PA® where we 

denote the system space with the subscript A and the 
ancilla with B. We are free to choose the initial ancilla 
state and so pick it to be w = |0) (0|, which gives as a 
solution 

{kB\U{t)\QB)=Mk{t) \jk. (2) 

This provides a dilation where the dimension of the an¬ 
cilla is R, and is known to be the smallest sized ancilla 
which may be required. The dilation unitary has dimen¬ 
sions Rd X Rd and Eq.([^ constrains d of its columns. 
These can be thought of as forming d orthonormal vec¬ 
tors in an Rd dimensional space. As we also need 
to ensure that the dilation is unitary, we impose that 
U{t)W{t) = I which is equivalent to requiring that the 
remaining d{R — 1) columns complete the orthonormal 
space. We prove in appendix A that this can always be 
done in a smooth way whenever the Mk(t) are smooth. 
The desired Hamiltonian satisfies the Schrodinger equa¬ 
tion, and is therefore given by H{t) = i U{t)W{t). 

The Hamiltonian involves the derivative of U, which 
is continuous but in general not analytic in t, thus there 
is the possibility of it being discontinuous or divergent 
due to the behaviour of the derivative of A 

careful analysis shows two potential problems. Firstly, 
Mkit) can be discontinuous (but always bounded) when 
Xkit 7 ^ 0) = 0. Secondly, it diverges when Afc(O) = 0 
but Afc(O) 7 ^ 0. These properties are inherited by the 
Hamiltonian, bar some accidental cancellation. The first 
case, where the Hamiltonian has a step change at some 
later time, corresponds to the Kraus rank of the chan¬ 
nel decreasing at a single point in time. More interest¬ 
ing is the second case. The divergence of the dilation 
Hamiltonian at t = 0, is avoided if and only if the dis¬ 
sipative part of e{t = 0) vanishes (see appendix B for 
proof). This gives the immediate corollary that all non¬ 
trivial time-independent Markovian channels lead to a 
divergent dilation. Such channels necessarily have that 
the survival probability of certain states decays linearly 
at short times. However, we know from the quantum 
Zeno effect that the survival probability of a state in any 
unitary system with a bounded Hamiltonian must de¬ 
cay quadratically for short times [30]. This linear decay 
might be a signature of the unbounded system environ¬ 
ment interaction, but often it is merely a consequence of 
approximations, like infinitely fast relaxation within the 
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environment na, made in the derivation of the Master 
equation. In both cases, the divergence of the dilation is 
inherited from this and is an indicator of Markovianity. 

Both the discontinuities and the divergence are, when 
they happen, benign. The divergence at t = 0 has only a 
finite impact on the dynamics as the dilation can always 
be picked such that U —)■ I as t —)■ 0. This implies that 
we can approximate the evolution arbitrarily well by re¬ 
placing the exact Hamiltonian by a bounded one. This 
can be seen from the relation between the error in the 
unitary (which is the error in the evolution) and that of 
the Hamiltonian [3T] : 

\\U{t)-UTim< t\\Hif)-HTit')\\dt\ (3) 

Jo 

where 11 • 11 is the operator norm and Ht and Ut are the 
target Hamiltonian and unitary respectively. We note 
that the right hand side is bounded by /p \\HT{t')\\dt' 
for any reasonable applied Hamiltonian, and is always fi¬ 
nite (even if Ht{0) is not). Thus the error in the unitary 
is always finite and can be decreased arbitrarily by hav¬ 
ing the applied Hamiltonian differ from the exact Hamil¬ 
tonian for a sufficiently short time. We calculate this 
error for a specific dilation later, and show that it can 
indeed be made arbitrarily small with ease. In the case 
of bounded discontinuities, the dilation Hamiltonian it¬ 
self can be arbitrary well approximated by a continuous 
one, which leads to the evolution being arbitrarily well 
approximated. This gives the following result: 


Hamiltonian. We start with the equation of motion in 
the given frame: 


f = (4) 

This can be represented in a different frame, for example 
in the interaction picture, by the transformation p(t) = 
Uo(typ{t)Uo{t) where Uo{t) = The equation of 

motion in this frame is given by: 

_ I = 

where: Lt(-) = U^LtiUo ■ U^)Uo + i[Ho, •]. 

We can relate the dilation of this Lindbladian, with uni¬ 
tary U(t)AB and Hamiltonian HAsit), to the dilation of 
the original Lindbladian. We first note that the respec¬ 
tive channels obey 


as p{0) = p(0). This implies that 


( 6 ) 


(7) 


which directly shows us that the dilation in the origi¬ 
nal frame is given by Uab = (?7o ® and so the 

Hamiltonian reads 


For an analytic family of quantum channels 
acting on d dimensional states there always 
exists a continuous and bounded Hamiltonian 
acting on at most d^ dimensional states such 
that the dynamics on the reduced system are 
arbitrarily well matched at all times. 


Further Methods 

The method we have just discussed works well, but it 
requires us to diagonalise the Choi state and, in the case 
that we are starting from an equation of motion, to calcu¬ 
late the channel. In practice, one or both of these may be 
very difficult to do analytically. Indeed, systems where 
these are hard to do are the ones we most want to find a 
dilation for and be able to simulate on a quantum com¬ 
puter, as they are precisely those which are difficult to 
simulate classically. We therefore present three alternate 
methods to construct dilations for complicated systems 
which rely on having found a dilation for a simpler sys¬ 
tem. The methods are: changing frames, separating into 
commuting parts, and perturbation theory. 

Firstly, we consider how the change of frame of the 
initial problem translates into a change in the dilation 


Hab=Ho®IF {Uo O I)Hab{U^o ® I)- (8) 


This gives us a very simple relation between a dilation 
Hamiltonian in one frame, and a dilation in a differ¬ 
ent frame. Indeed, the relation is the same as when we 
change frame in normal unitary dynamics with the small 
additional step of transforming Uq —>■ Uq ^ I. 

Secondly, we show how a dilations can be calculated by 
separating the Lindbladian into commuting parts, at the 
cost of increasing the ancilla space. Take two different 
channels which commute at all times such that 

4">(0 = <!l‘’(To). (9) 

By performing the dilations one after the other, and on 
different ancilla spaces, we have that 




^ ab 


) 


( 10 ) 


= Trsc 


rr(l)rr(2) (X) 

'^AC^AB PA‘i 9 UJB ^AB^AC 


where Uac = Uac ® Is and similarly for the other op¬ 
erators. This gives as the Hamiltonian 


77(12) _ prfi) I rr(i) 

^ABC ~ ^ An A- u A nOi ^ 


^AC 


^AC^^AB'-'AC ■ 


( 11 ) 
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Due to the commutativity, we can inverse the order in 
Eq.@ and obtain 

^ABC — ^AB + ^AB^AC^AB 

which is, in general, a different Hamiltonian but leads to 
the same dynamics on the reduced system. 

Thirdly, we can use a perturbative approach. We take 
the Lindbladian for the system to be with 

i5 <C 1, where we assume we have already found a dilation 
for and find a new dilation which gives the correct 
dynamics to first order in <5. The quantum channel for 
such a Lindbladian is 


Examples 

We demonstrate these methods and results by look¬ 
ing at specific examples of quantum channel acting on 
qubits which represent some of the iconic decoherence 
models in quantum information [UlllllMEl- In many 
cases the resulting Hamiltonians are sufficiently simple to 
be directly experimentally realisable. Consider the spin- 
boson model, where a single spin interacts with a bath of 
bosons via a constant Hamiltonian, such that the reduced 
dynamics are given by p = [cr^jp]], where the 

decay rate is a function of time [351ISS] ■ Performing the 
steps outlined above, we find the dilation to be 


et — f-t 
^ .(0) 


^(0) ;-(l),(0) 




0{S^) 


Se. 


( 1 ) 


(13) 


where Constructing the Choi state 

from the channel is a linear transformation which can be 
done separately for and The eigenvalues and 
vectors of -|- can be found to first order in <5 us¬ 
ing standard methods from perturbation theory. This is 
much easier to do than diagnosing the Choi state exactly, 
although some of the advantage is lost if the original state 
had a high degree of degeneracy which is broken. Hav¬ 
ing done this, we can easily find the Kraus operators by 
expanding Eq.(l) in the main body to first order in S. It 
is worth noting that if = 0 it is sufficient to find 

Mk to O(v^), as the equation of motion is quadratic in 
the Kraus operators and we are only interested in find¬ 
ing the dynamics up to first order in S. In the case that 
the Kraus rank of the channel is unaffected by the per¬ 
turbation, so that the above condition does not hold, we 
can write Mk = In that case the 

correction to the dilation unitary satisfies 

{kB\U^^Ht)\0B) =Mi^\t) Vfc. 

[/(i){7(o)t + [7(o)[/(i)t = 0. (14) 


The last expression is equivalent to requiring the unitar- 
ity condition to hold to first order in <5. This reduces 
the problem of finding the dilation unitary to a system 
of linear equations. The Hamiltonian is then given by 




such that H = -|- + 0{S^). In the case that 

the Kraus rank of the channel does change, similar ex¬ 
pressions can be found, although care must be taken to 
ensure that both terms 0{\/S) and 0{d) are properly ac¬ 
counted for. This perturbative method can be extended 
to take into account second, or higher, order effects. 


Hs.bit) 


lit) 




(16) 


where CTj, acts on the system and ay acts on the ancilla 
(which consists of a single qubit). The exact form of the 
Hamiltonian depends on 7(1); for a typical spectral den¬ 
sity we have that 7(0) = 0 and that, for certain values 
of t, it becomes negative. In such a case the channel is 
non-Markovian and the dilation is bounded and continu¬ 
ous for all times. In the case of constant 7, however, this 
is a Markovian dephasing Lindbladian which gives rise 
to the phase-flip channel, and whose dilation diverges at 
t = 0 as expected. If we approximate this Hamiltonian 
to be constant by replacing the scalar pre-factor by C 
for short times (that is, when > tJ) the error 

between the target unitary and the unitary reached is 

2 

upper bounded by ^ + O(^). Thus, with a sufficiently 
large C the error can be made arbitrarily small. Another 
way of looking at the errors is to see how the dynam¬ 
ics of a state depends on C, as is plotted in FIG.(|^. 
We see that the main effect of introducing a cut-off is 
to make the behaviour quadratic, rather than linear, at 
short time, and that even a modest value for the cut-off 
is enough to reduce errors to insignificance. 



■ Cut off =5y 
-- Cut off =25}' 
— Exact 


FIG. 2. The decay of the |-|-) (-|-| state when subjected to a 
dilation of the dephasing channel, where the divergent Hamil¬ 
tonian is replaced by a finite cut-off at short times. The inset 
shows the behaviour for short times. 

As another example, we look at the dilation for an 
amplitude damped system |3D]. Here a non-degenerate, 
two-level system relaxes to its ground state in a process 
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such as spontaneous emission. We start with an equation 
of motion p = — 7 ({(T+cr_,p} — 2 cr_pcr+) + where 

a± is the raising/lowering operator, and find the dilation 
to be 


Ha.d.{t) = 


I'y 


- 1 


(cr_ ® — a+ ® a-) + ® I, 


(17) 

where once again the ancilla consists of a single qubit. 
It is interesting to note that, although this has a very 
different physical origin to the dephasing channel, the 
time dependency is almost identical. The comparison 
highlights some of the most common features of dilations 
of simple systems. The Hamiltonian is strongest at t = 0 
(possibly even diverging, as discussed previously), and 
the terms corresponding to decay fall to 0 for large t; 
which is to be expected for the system to settle to its 
steady state. 

As a more involved example, we now consider a qubit 
acted on by an amplitude damping Lindbladian as above 
but with an additional constant driving term, —iU,[ax, p]- 
In this case the coherent and incoherent part of the evolu¬ 
tion no longer commute, which makes a direct calculation 
of the exact dynamics difficult. Nevertheless, by taking 
the limit where the driving strength H is much smaller 
than the decay rate 7 , we can use perturbative methods 
to find the dilation 


H{t)=Ha.d.{t)+n 


1 -beT'* 


CTa; 0 I -b H 


y/ — I 
'(eT-* -b I)^^ 


( 18 ) 

where we have set Wg = 0 for simplicity. This dilation has 
two potentially unexpected features. Firstly, the driving 
term has gone from being constant to decaying in time. 
For large t it does this at same rate as Ha.d., which is 
necessary for the system to settle to a single fixed point. 
Secondly, we have the appearances of a third term, which 
is back action caused by dilating the control and it has a 
more complex structure in time, although it too decays at 
the same rate for large t. This term is caused by Ha.d. {t) 
building up entanglement between the system and the 
ancilla. The emergence of complicated time structure in¬ 
duced by a simple control field is closely linked to the fact 
that Master equations are changed in a non-trivial way 
by the addition of a Hamiltonian acting on one subsystem 


m- 


A more complex case involving instead a time depen¬ 
dent driving term, —ifl cos(u;t) [cTj,, p], can also be dilated 
in a perturbative method. In this case we also make the 
Rotating Wave Approximation and, in the resonant case, 
the dilation is: 


H{t) = iHo{t) a-^ a++ (19) 

-b ft, f{t) cr_ 0 1 -b g{t) Uz (Si (Jx + h.c. 


where Ho{t) = e — I, f{t) = 

e“*““‘/(I + e"''*), and g{t) = — 1/(4 (e'*'* -b I)^) are 



FIG. 3. The dilation for a qubit subjected to an amplitude 
damping Lindbladian and a resonant sinusoidal driving field 
is detailed in Eq.(19l. We plot here the real part of the 
time-dependent functions in that Hamiltonian with = 2. 
Only one of the terms diverges at t = 0, while all the terms 
decay to 0 quickly at large times. 


plotted in FIG§ The increase in complexity of the 
Hamiltonian is directly related to the time-dependence 
of the original equation of motion, but the dilation could 
still be constructed which shows that this approach is 
applicable to a wide range of problems. 


Conclusion 

By rescaling time in a nonlinear way, some of these 
channels can even be dilated to constant Hamiltonians. 
In general, this is true whenever the dilation Hamiltonian 
is of the form Hit) = h{t)X such that it has a single time 
dependent factor. This allows us to apply a constant 
Hamiltonian H = h^X for time r to simulate the real 
dynamics evolving for a time t where t Jq h(t')dt'. 
As h{t) is continuous and bounded, this is always well 
defined. In cases where there are several different time 
dependencies, this method can be used to remove one 
of them. This is particularly useful in eliminating diver¬ 
gences, which would otherwise be problematic to imple¬ 
ment experimentally. It also has the advantage that, in 
many cases, the evolution for an infinite amount of time 
t can be simulated with a finite r. 

As well as being useful in the single qubit case, these 
results can also be directly used in the case of a sys¬ 
tem of many qubits all subjected to independent chan¬ 
nels, leading to a superpolynomial speed up from what 
could be achieved classically. Prominent questions which 
could be investigated include the decay of multipartite 
entanglement of an initially highly entangled state |38j . 
or the performance of one-way computation [39) with a 
cluster state affected by local noise. For such systems 
of N qubits the dilation can be calculated once classi¬ 
cally and the dynamics of the entire system then simu¬ 
lated on a quantum computer using a maximum of 2N 
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ancilla qubits and Hamiltonians that affect, at most, 3 
qubits. Simulating this classically would require apply¬ 
ing the quantum channel up to 2^ times for an initial 
state which is highly entangled or, equivalently, solving 
the dynamics of the complete channel which would have 
2^ Kraus operators. 

This method also provides the ability to do something 
which a normal Stinespring dilation cannot do at all. By 
simulating the dynamics continuously in time, the evolu¬ 
tion of the state follows the ‘true’ path that we are trying 
to simulate, rather than just reaching the required goal 
for a single instance in time. This means that the full 
information about the behaviour of the system over an 
extended interval of time is accessible, allowing simula¬ 
tions where the time at which measurement takes place 
is not known a priori. In our scheme, such a scenario 
creates no additional difficulty, as the system follows the 
correct continuous dynamics. In a standard Stinespring 
dilation approach, however, this is either completely im¬ 
possible or introduces substantial errors. These can be 
quantified by considering the snapshot Stinespring dila¬ 
tion as being a constant Hamiltonian (the logarithm of 
the unitary dilation) applied for a shorter or longer pe¬ 
riod of time, leading to an error which grows in time. 

In our approach, the complexity of open system dy¬ 
namics is condensed in the time-dependent system envi¬ 
ronment interaction, allowing a simulation to be imple¬ 
mented using state of the art methods; whereas in nat¬ 
ural systems, and in previous approaches to simulating 
the dynamics of open systems, the complexity resides in 
the dynamics of the infinitely large environment. These 
two perspectives can be understood as the two ends of 
the spectrum of quantum simulations. Since any time- 
dependence can be understood as originating from the 
dynamics of an additional system (fundamentally all in¬ 
teractions are time-independent), our approach suggests 
very clearly how to access the entire spectrum: expand 
the ancilla system while gradually reducing the time- 
dependence of the interaction. Such a continuous varia¬ 
tion will give valuable insight, for example in the control¬ 
lability of open quantum dynamics, as seen in Eqs.(18) 
and (19). There the back action of an external con¬ 
trol caused by the environment interaction makes itself 
transparent, whereas such effects are extremely hard to 
unravel in a model based on an infinitely large environ¬ 
ment. The suggested transition would allow the study 
of this back-action in its entire range of manifestations, 
opening up a completely new angle on the investigation 
of control on open quantum systems. This, in turn, fos¬ 
ters our endeavours in the struggle against decoherence 
and the realisation of working quantum technologies. 
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Appendix A: Completing the Unitary 


due to the orthogonality at t. This explicitly shows that 
|fd+i(^)) is a smooth function of t. From this we can 
easily obtain the normalised vector 

\vd+i{t + 5t)) = \vd+iit)) (23) 

n—d 

+ 6t'^Re {{vd+i{t)\An{t)) {vd+i{t)\vnit + 6t))) \vd+i{t)) 

n—1 

n—d 

{vd+l{t)\An{t)) \Vn{t + St)) + 0{5i^) 

n—1 

which is also a smooth function of t whenever the Mk{t) 
are smooth. By induction, we see that this is true for 
all \v{t)) and all t. Hence we can always construct a 
unitary which is smooth whenever the Kraus operators 
are smooth. 


In order to find a Hamiltonian for a dilation, we first 
have to find the corresponding unitary. As noted in 
Eq.(2), this unitary satisfies 

{kB\Uit)\0B)=Mkit) Wk, (20) 

C/(f)C7l'(t) =1. 


We show that it is always possible to find a solution for 
U{t) which is smooth whenever the Mk{t) are smooth. 
The first condition constrains d columns of the unitary 
as being orthonormal vectors and the second condition 
requires us to pick the remaining d{R — 1) columns such 
that they form a complete orthonormal set. At t = 0 we 
use Gram-Schmidt from an arbitrary basis such that the 
none of the d{R— 1) columns are 0. In order to ensure at 
later times that these extra columns vary continuously in 
time we use the vectors at t as a basis for Gram-Schmidt 
at t + St. Let the columns of U{t) be expressed as the Rd 
vectors \vn{t)). At time t + St the first d of these vectors 
are transformed according to the change of the Kraus 
operators. They are still orthonormal to each other, but 
no longer orthogonal to the other R{d — 1) vectors. Us¬ 
ing Gram-Schmidt, we update the d -I- 1 vector by first 
calculating: 

\vd+i{t + = (21) 

n—d 

\vd+l{t)) - ^ {vd+l{,t)\Vn{t + St)) \Vn{t + St)) . 

n—1 

For n < d we can also write 


\Vn{t + St)) = \Vn{t)) + St |A„(t)) -I- 0{Sf) (22) 


whenever the |u„(f)) are a smooth function of t, that 
is, whenever the Mkit) are smooth. This enables us to 
rewrite Eq.(21) as 


\v’d+^{t + St)) = 

n—d 

\vd+i{t)) -'^St {vd+i{t)\An{t)) |■U„(^ -k St)) + 0{St^) 

n—1 


Appendix B: Divergence at t = 0 

We prove that the dilation of e(<) diverges at f = 0 if 
and only if e(0) is non-Hermitian. We first consider the 
case where the dilation does not diverge, which lets us 
write 


e* (p) = Ttb [U{t)p (g) wUl (t)] 

et{p) = Trs {-i[H{t),Uit)p(S)UjUHt)]} (24) 

eo(/o) = -iTrs {[iJ(0), p (g) w]} 

as eo(p) = p means that there exists a dilation such that 
U (0) = I. We expand H (0) in terms of separable Hermi- 
tian operators: 


H{0) = Y,Ak^Bk. (25) 

k 

This enables us to calculate the partial trace according 
to 


eo(p) = -*TrB 


{Akp (g BkUj - pAk (g ujBk) 


= -»y^[AfcAfc,p] 


(26) 


= -i[H',p\ 


where Afc = Tr^uiBk). These are necessarily positive as 
Bp is Hermitian and w is a state, from which it follows 
that H' must also be Hermitian. This proves that if the 
dilation is finite at t = 0, then eo(p) must be purely 
Hermitian. 

To prove the converse we recall that the divergence 
depends on the eigenvalues of the Choi state. Specifically, 
that at least one of the eigenvalues must obey A(0) = 0 
but A(0) yf 0. The Choi state is pure if and only if 
the quantum channel is unitary. If this is the case it is 





clear that the eigenvalues do not change, so we ignore the 
Hermitian part of the channel. Assuming that cq = 0, we 
can write the quantum channel for short times as 

et«I + A^+0(t3). (27) 

As previously stated, reshuffling the elements of Ct, writ¬ 
ten as a matrix, gives us the Choi state A(. This means 
that the elements of At have, potentially, terms of every 
order in t apart from first order. As the Af is a state, 
though not quite normalised, its eigenvalues are all non¬ 
negative and sum to d. It is also a pure state at t = 0 
and, as the eigenvalues are analytic in t, we can we write: 

Aq = d -|- cxot -\- 0{t^) 

Xi = 0 + ait + + 0{t^) i = — l (28) 

The eigenvalues are also obtained from the characteristic 


equation for At, and so we can write: 

Xd^+AX‘^^-^ + BX‘‘"-^ + ... = 0 (29) 

where A, B, C,... are products and sums of the elements 
of At. This means that none of these coefficients can have 
terms which are linear in t. From the basic properties of 
the roots of polynomial equations we have that 

B = '^X^X^ pL,v = Q,...,d? - I 

= d'^^ait + 0{t‘^). (30) 

i 

But, as we have already noted, B cannot have a term 
proportional to t. This implies that oti = 0 and, as 
each are non-negative, that at = 0. The sum of the 
eigenvalues hence also means that ao = 0- Therefore, 
all the eigenvalues obey X{t = 0) = 0 and, as noted 
previously, this is a sufhcient condition for the dilation 
to not diverge. 



